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The symmetry operations of the crystal groups relevant for the high temperature superconductors 
HgBa 2 Cu0 4 +« (Hgl201), YBa 2 Cu 3 07-<s (YBCO), and BiaSraCaCuaOs+a (BSCCO) are elucidated. 
The allowable combinations of the superconducting order parameter (OP) components are presented 
for both the angular momentum and lattice representations. For tetragonal Hgl201, the spin singlet 
OP components are composed from four sets of compatible basis functions, which combine to give 
the generalized s-, d x 2_ y 2-, d xy -, and g xy ( x 2_ y 2)-wa,ve OPs. Although YBCO and BSCCO are 
both orthorhombic, they allow different mixing of these components. In YBCO, the elements of 
those s- and d x 2_ y 2-wa,ve sets (and of those d xy - and g xy ( x 2_ y 2ywwe sets) are compatible, but in 
BSCCO, elements of s- and d xy - wave sets (and of d x 2_ y 2- and g xv ( w 2_ v 2\-wsve sets) are compatible. 
The Josephson critical current density j/ across c-axis twist junctions in the vicinity of T c is then 
evaluated as a function of the twist angle cj>o, for each allowable OP combination, for both coherent 
and incoherent tunneling. For tight-binding Fermi surfaces, we argue that coherent tunneling is 
only possible for cj>o = 0°,90°, for which the rotated Fermi surfaces line up. Recent experiments 
of Li et al. demonstrated the independence of j£ ((j>o) / Jc upon <f>o at and below T c , where J,f 
is the critical current density of a constituent single crystal. These experiments are shown to be 
consistent with an OP containing an s-wave component, but inconsistent with an OP containing the 
purported d x 2_ y 2-~w&ve component. In addition, we argue that they demonstrate that the interlayer 
tunneling across untwisted layers in single crystal BSCCO is entirely incoherent, with no amount of 
forward scattering. We propose a new type of tricrystal experiment using single crystal c-axis twist 
junctions, that does not employ substrate grain boundaries. 
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I. INTRODUCTION 

There have been some exciting new developments 
which reopen the question of the symmetry of the order 
parameter (OP) in the high temperature superconduct- 
ing compounds (HTSC). Although previous exper- 
iments such as the penetration depth, and "phase- 
sensitive" tricrystal experiments were interpreted as giv- 
ing strong evidence for an OP in YBa2Cu307_«5 (YBCO) 
that was dominated at low temperature T by a compo- 
nent of d x 2_ y 2 symmetry, other experiments showed 
that there was also a substantial s-wave component in 
that material as well, |)],[l0| and a very recent experi- 
ment suggested a nodeless OP. || However, as discussed 
in the following, the crystal symmetry of YBCO is or- 
thorhombic, with distinct Cu-0 bond lengths a and b 
within the Cu02 planes, which are normal to and along 
the CuO chain direction, respectively. This particular 
type of orthorhombicity then allows for an arbitrary mix- 
ing of d x 2_ y 2- and s-wave OP components, even at T c . 

To date, the only hole-doped cuprate superconduct- 
ing material known to have a consistently tetragonal 
crystal structure is HgBa 2 Cu04 + 5 (Hgl201). Unfortu- 
nately, no "phase-sensitive" experiments have yet been 
performed on this material. The only experiment rele- 
vant to the symmetry of the OP was point contact tun- 
neling, which appeared consistent with an s-wave OP at 



low T. |Tl]] A material which can sometimes be made in 
a tetragonal structure is Tl 2 Ba2Cu04 +( 5 (T12201). Point 
contact tunneling and Raman scattering on T12201 ap- 
peared to be consistent with line nodes of the OP, as 
expected for a d x 2_ y 2-w&ve OP. |jl^,[l3) There was one 
phase-sensitive experiment performed on T12201, which 
claimed to give strong evidence for a d x 2_ y 2--w&ve OP in 
that material at low T. (m). However, T12201 is diffi- 
cult to prepare in the tetragonal form, and most sam- 
ples of it are actually orthorhombic, in a form simi- 
lar to that of YBCO. jl5| The electron-doped cuprate 
Nd2- rc Ce a; Cu04_5 (NCCO) is also tetragonal, and mi- 
crowave experiments suggested a full BCS-like gap, |i"6| ] 
suggestive of conventional s-wave superconductivity in 
that material. However, new Pb/NCCO c-axis Josephson 
tunneling experiments demonstrate a surprising small s- 
wave OP component, and no evidence of a gap, suggest- 
ing that the full OP might not be pure s-wave. || 

With regards to Bi 2 Sr 2 CaCu 2 8+ 5, (BSCCO), recent 
phase-sensitive tricrystal experiments also provided ev- 
idence at the tricrystal intersection point for a half- 
integral flux quantum $0/2, where the standard su- 
perconducting quantum of flux is $0 = hc/2e. filf | 
In addition, angle-resolved photoemission spectroscopy 
(ARPES), penetration depth measurements, point con- 
tact tunneling measurements, and Raman scattering ex- 
periments were all interpreted as providing additional ev- 
idence for a superconducting OP consistent with d x 2_ y 2 
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symmetry, although such non-phase-sensitive experi- 
ments could only detect the magnitude of the OP. 1 18 2^] 

Early attempts to observe c-axis Josephson tunneling 
between Pb and BSCCO were unsuccessful. p3fl . Those 
experimenters prepared thin films of BSCCO, and de- 
posited a thick layer (roughly 1000A) of Ag on top of it, 
followed by a thicker layer of Pb. However, new c-axis 
Josephson tunneling between Pb and BSCCO provided 
compelling evidence for a small s-wave OP component at 
low T. Q These experimenters cleaved a single crystal of 
BSCCO, and deposited only about 10 A of Ag on it prior 
to depositing Pb. The product of the c-axis supercon- 
ducting critical current I c and the quasiparticle resistance 
R n was found to be very small, about 1-2 /iV, curiously 
about the same as for the new Pb/NCCO junctions. || 
Nevertheless, this evidence for (at least) a small s-wave 
OP component in BSCCO at low T is very interesting, 
since BSCCO is orthorhombic in a different sense than 
is YBCO. As discussed in the following, this form of or- 
thorhombicity does not allow for a mixing of d x 2_ y 2- and 
s-wave symmetry at T c . Thus, these c-axis Pb/BSCCO 
Josephson experiments could only be explained in the 
"d-wave scenario" (in which the purported d x 2_ y 2-w&ve 
OP component were dominant at T c ) by the appearance 
of a second phase transition in BSCCO between the T c 
values of BSCCO and Pb. The same argument would of 
course apply to the ei-wave scenario as an explanation of 
the new Pb/NCCO experiments. |6| 

Furthermore, c-axis twist experiments have been per- 
formed, in which a single crystal of BSCCO was cleaved 
in the a6-plane, the two single crystal pieces twisted an 
angle 4>o with respect to each other, and fused together 
just below the melting temperature. Preliminary ex- 
periments performed at low T and large magnetic fields 
H||c showed that the measured I c across the twist junc- 
tion did not depend upon 4>q. However, these early re- 
sults could be interpreted in terms of a dominant d x 2_ y 2- 
wave OP component that could "twist" by mixing in a 
subdominant d^-wave OP component. pi| In the vicin- 
ity of a 45° twist junction, the dominant d x 2_ y 2--wave OP 
component would be suppressed, allowing the subdomi- 
nant d xy -wave OP component to become non-vanishing 
well above T C 2, even extending to some extent up to T c . 
Thus, one could have a rather (^-independent / c (</>o) be- 
havior in the vicinity of the measurement T = 10 K, as 
observed. [Q However, such a scenario would necessarily 
also imply a second bulk phase transition at T C 2 <<T C , 
below which the e^y-wave form could appear. T C 2 would 
have to be very low (i. e., T C 2 < IK) in the bulk of 
the sample, if the states in the superconducting "gap" 
evident at low T were correctly interpreted as arising 
from nodes in the purported d x 2_ y 2-wa.ve OP. |l8|-p2|] In 
addition, although such a scenario could lead to a non- 
vanishing / c (0o = 45°) in the vicinity of T c , it could not 
lead to an isotropic I c (<po) near to T c . 

Now the c-axis twist experiments have been performed 



just below T c in BSCCO. |2|J3|] These new experiments 
also show no 4>q dependence of I c ((f>o). In these exper- 
iments, I c was found to scale with the twist junction 
area, indicating uniform junctions, and the critical cur- 
rent density J/ of the twist junctions was found to be 
the same as , the critical current density for the con- 
stituent single crystals, regardless of 0o- Since BSCCO 
is intrinsically a stack of Josephson junctions, [^5| these 
experiments demonstrate that the twist junctions behave 
precisely as the Josephson junctions intrinsic to single 
crystal BSCCO. Based upon our earlier detailed calcu- 
lations, these experiments are incompatible with the "d- 
wave scenario." [ p4[ Here, we present general arguments, 
based solely upon the crystal symmetry, which demon- 
strate the basic incompatibility of a d x 2 _ y 2 -wave with an 
s-wave OP component. As shown in the following, the 
experiments of Li et al. provide strong evidence that 
the OP in BSCCO contains the s-wave component at 
T c . We employ group theory to argue that this experi- 
ment demonstrates that the dominant OP at T c contains 
the s-wave component, but not the purported d x 2_ y 2- 
wave component. These new phase-sensitive experiments 
are thus in direct contradiction with the phase-sensitive 
tricrystal experiments of Tsuei and Kirtley et al., which 
were claimed to provide strong evidence for a dominant 
d x 2_ y 2-wave OP in BSCCO at low T. @ 



II. CRYSTAL SYMMETRIES 

Most workers in the field of HTSC agree that the su- 
perconducting properties of the materials are determined 
by the structure of the Cu02 planes. This is because the 
electronic structure calculations indicate that the Cu02 
plane bands cross the Fermi energy Ep , [^6| and ARPES 
measurements well above T c provide strong support for 
these results. [g7| Within the Cu02 planes, the Cu and 
O ions are situated approximately as pictured in Fig. la. 
p8| The Cu ions are located on a square lattice, with O 
ions between near-neighbor Cu ions. 

If one ignores the complications of c-axis variations in 
the relative orientations of different Cu02 planes, the 
crystal point group describing a single, tetragonal, Cu02 
plane is C/± v . [p9|-|33f The C^ v group symmetry operations 
are indicated in Fig. lb. In addition to the identity 
operation E, there are mirror reflections o~ x and a y in 
the planes normal to the layers and containing the x and 
y axes, respectively, mirror reflections a^i and 0^2 in the 
planes normal to the layers and containing the diagonals 
d\ and di, respectively, and rotations C4, C7 1 , and C2 = 
Cf by 90°, -90°, and 180° about the c-axis, respectively. 
Tinkham used the slightly different notation a v , o~' v , ad, 
and a' d , for the four mirror reflections, respectively. p9| ] 
This crystal symmetry is appropriate for Hgl201, NCCO, 
and for tetragonal samples of T12201. 32 
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In YBCO, there is an orthorhombic distortion involv- 
ing the CuO chains, which causes the CuC>2 crystal struc- 
ture to be longer along the b- (or y-) axis than along the 
a- (or x-) axis. The crystal space group is Pmmm, with 
the point group of the CuC>2 planes being C\ v , an exam- 
ple of Ci v . For this lower symmetry, the allowable group 
operations in addition to E are only a x , a y , and C2, as 
sketched in Fig. lc. This group is also appropriate for the 
double chain compound YE^C^Og (Y124). However, 
there are some interesting distortions, which complicate 
this structure. These are the buckling of the CuC>2 planes 
along the c-axis direction. ]3l| Nevertheless, there isn't 
any buckling within the CuC>2 planes, so it doesn't affect 
the symmetry arguments relevant to the CuC>2 planes. 

In order to correctly interpret the new experiments on 
BSCCO, it is imperative that the crystal symmetry be 
understood as thoroughly as possible. The discussion of 
the symmetry of the superconducting OP in this paper 
depends crucially upon the notion that there is at least 
one mirror plane that intersects the CuC>2 plane. This 
mirror plane then distinguishes OP basis functions that 
are odd or even with respect to reflections in it. As dis- 
cussed in the following, the 6c-plane is a strict crystallo- 
graphic mirror plane in most samples. 

For overdoped or optimally doped BSCCO, the crys- 
tal structure is also orthorhombic, but in this case the 
inequivalent a- and 6-axes are not along the Cu-0 bond 
directions, but along the directions di and d2 diagonal 
to them. In addition, there is a periodic lattice distor- 
tion Q along one of these diagonals, which we identify as 
the b- (or d 2 _ ) axis, consistent with the identification of 
most workers, [B4 35 although some workers identified 



Q as being along the a- (or d\-) axis direction. [ p6|]37| 
More precisely, Q is incommensurate along the recip- 
rocal lattice wavevector b*, but commensurate in the 
other crystal directions, Q = (0,0.212,1). The 
fact that there is a non-vanishing commensurate c-axis 
component actually removes the <id\ mirror plane crys- 
tal symmetry, consistent with the appearance of weak 
(100) reflections in the x-ray diffraction data. p5| , ^6| . 
In addition, there is some dispute as to the precise de- 
tails of the periodic lattice distortion. An early electron 
diffraction study of BSCCO indicated that in one sample, 
Q = (0.004,0.212,1), containing a small a* component 
that would also remove the &d2 mirror plane symmetry, 
p4| although this slight a* component of Q was not seen 
in high resolution x-ray diffraction experiments on re- 
lated samples. |34|] More to the point, electron diffraction 
studies of the samples used in the present c-axis twist ex- 
periments demonstrated that the small a* component of 
Q was absent in 90% of the samples. [B8j Thus, it appears 
that although the crystal symmetry of BSCCO is lower 
than orthorhombic, as the mirror plane Od\ is not pre- 
cise, most presently-made BSCCO samples exhibit the 
<jd2 mirror plane (the 6c-plane). Further experiments to 
clarify this situation are under way. |3§[] 



This periodic lattice distortion also gives rise to some 
buckling of the lattice, including Cu-0 bond buckling in 
the &c-plane. Both x-ray studies of the crystal struc- 
ture of BSCCO used a four-dimensional analysis proce- 
dure, and listed the basic space group as Bb2b (or equiva- 
lently, A2aa) , which is the orthorhombic point group , 
another form of C^v Unfortunately, they differ 

slightly in listing the full crystal group as M : A2aa : 111 
and Niii h , respectively. In addition, there is a fur- 
ther complication in some heavily underdoped samples, 
in that they appear to undergo a phase transition from 
orthorhombic to monoclinic. |3^,Q In such cases, there 
are no mirror planes (other than parallel to the a6-plane) 
remaining, and the only non-trivial group operation ap- 
propriate for a single Cu0 2 plane would be C 2 . How- 
ever, since the c-axis twist experiments did not involve 
such monoclinic single crystals, Jl|,^8| we shall not discuss 
such very low symmetry cases further. 

Neglecting these fine details, the allowed symmetry op- 
erations of the simplified point group are E, o~di, 
(Jd2 1 and C2, as illustrated in Fig. Id. This symmetry 
is also appropriate for the compound L^-zSr^CuO^a 
(LSCO). With regard to a single Cu02 layer, the Od\ 
mirror-plane symmetry is strict, as it is only broken by 
the c-axis component of Q in BSCCO. More important, 
C\ v (for YBCO) and Cf v (for BSCCO) differ only in that 
the in-plane crystal axes are effectively rotated by 45° 
about the c-axis from each other. As regards the orbital 
symmetry of the superconducting OP, this difference is 
crucial as noted in detail below. 



III. ORDER PARAMETER COMPONENTS 

The role of symmetry upon the OP components in the 
HTSC was discussed in a recent review article by Annett, 
Goldenfcld, and Leggett (AGL). (3^] However, the nota- 
tion used by AGL was rather confusing, such as using 
"s~" to denote a state that contains no s-wave compo- 
nent, and whose leading term is g- wave. More important, 
the question of whether different OP components could 
mix with or without a second phase transition was stated 
incorrectly in the conclusion. Therefore, as the symmetry 
of the OP in the HTSC is still very much under debate, 
it is imperative to provide a clear discussion of the issues. 

Quantum mechanics, Bloch's theorem, and group the- 
ory dictate that the superconducting OP must be consis- 

Assuming that the 



tent with the crystal structure. |29 



Cu02 planes are responsible for the superconductivity in 
the HTSC, then the point group structure of the Cu02 
planes will identify those OP components that are com- 
patible or incompatible with each other. For a tetragonal 
crystal of Hgl201 or NCCO, these OP components act 
under the point group C^ v according to Table I. The OP 
components fall into sets, each of which has the same set 
of eigenvalues in the table. Each set of OP components 
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with the same eigenvalues under the group operations 
forms a basis, from which an OP eigenfunction can be 
constructed. We have labeled these OP eigenfunctions 

as |s), \d x 2_y2), \d X y) , \9xy(x 2 —y 2 ) 

), and \ (p x ,p y )), respec- 
tively. We have also listed the standard group theoretic 
(GT) notation for each OP basis set. [^9| The spin singlet 
OPs have one-dimensional representations, with eigen- 
value +1 under E. Those labeled A and B are even and 
odd under C4, respectively, and those with subscripts 1 
and 2 are even and odd under a x and a y . Spin triplet 
OPs have two-dimensional representations, with eigen- 
value +2 under the identity operation E, and thus have 
GT notation E. H|l£l) 

But, how can we best represent the OP basis sets 
mathematically? We shall employ two methods. In 
the usual BCS theory, pairing occurs between otherwise 
free quasiparticlcs with opposite momenta ±k^ on the 
Fermi surface. For a free-particle Fermi surface, the mo- 
menta are then described by the planar angle 0k, where 
ki? = kp (cos 0k, sin 0k). The OP basis is then described 
as in scattering theory, in which the scattered waves are 
expanded in terms of angular momentum quantum num- 
bers £. On the other hand, for real-space pairing, the 
pairing takes place between quasiparticles on distinct lat- 
tice sites, such as on the same or nearest-neighbor sites. 
One then writes k = (k x ,k y ), which components may 
satisfy a relation on the real-space Fermi surface. The 
first Brillouin zone (BZ) is given by \k x \, \k y I < ir/a. The 
real-space quantum numbers, count the lattice 

separations between the paired quasiparticles along the 
Cu-0 bond directions. 

In either scenario there are an infinite number of pos- 
sible OP basis elements, since I is any integer satisfying 
I > and n, m can be any integers. However, the crystal 
symmetry forces the infinite number of OP basis elements 
to appear together in a limited, finite number of OP basis 
sets, each set having an infinite number of elements. We 
thus describe each possible OP as having a set of com- 
ponents, each component represented by a basis element. 
Each element in the OP basis set transforms identically 
under each of the allowed crystal point group operations. 

In the angular momentum £ representation, each OP 
basis function (component) may be written as either 
cos(^0k) or sin(^0k)- In Fig. 2, we sketch the even £ 
OP basis functions for I < 4 appropriate for singlet spin 
pairing. Shaded and unshaded regions describe oppo- 
site signs of the function. Shown are representations of 
the I — s-wave function, a constant, the two £ — 2 
functions d x 2_ y 2 and d xy , proportional to cos(20k) and 
sin(20k)> respectively, and the two I = 4 functions g x 2 y 2 
and g xy ( x 2- y 2), proportional to cos(40k) and sin(40k), 
respectively. 
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TABLE I. Superconducting OP character table for a 
tetragonal crystal such as Hgl201, with point group Cav- The 
symmetry operations are the identity (E), reflections about 
the x (<T X ), y (o~y), di (a d i), (o"d2) axes, and rotations about 
the c-axis by ±90° (C 4 , C7 1 ) and 180° (C 2 ). The group theo- 
retic (GT) notations Ai, A2, B\, and B2 are one- dimensional 
representations, and E is a two-dimensional representation. 

The symmetry operations of group Ca v (Table I) di- 
vide the even £ OPs into four basis sets, which are listed 
in Table II. These sets consist of the basis functions 
cos(^0k) and sin(£0k), where £ is even, and either is or 
is not divisible by 4, as indicated in Table II. Thus, the 
g x 2 y 2 function, cos(40k) (equivalent to 1 — 8fc^/c^, which 
averages to zero), is an element of the same OP ba- 
sis set (|s)) as the £ — 0, s-wave function, a constant. 
It is possible to define general s-wave eigenfunctions, 
Pi(0k) = S^Lo -^n cos(4n0k), each of which transforms 
as the basis set denoted A\ in Table I. We are then able 
to rewrite each of the OP eigenfunctions in terms of the 
lowest £ basis functions multiplying one of these general- 
ized s-wave eigenfunctions, as shown explicitly in Table 
II. 



OP Eigenfunction £ Representation 
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TABLE II. Angular momentum quantum number £ 
representations of the singlet superconducting OP eigen- 
functions for a tetragonal crystal. The functions 

r *(<M = E:io^ cos ( 4n< M- 
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In the lattice representation, the general tetragonal 
eigenfunctions are cos[a(nk x + mk y )] and sva[a(nk x + 
mky)], where n and m can be any integers. One clas- 
sifies the OP basis functions in terms of the spatial sep- 
arations of the paired quasiparticles. In order of in- 
creasing separation, there is on-site (for a retarded BCS- 
type interaction), near-neighbor, next-nearest-neighbor, 
etc. pairing. Such expansions of the OP eigenfunctions 
are presented in Table III. In this representation, the 
most general OP eigenfunctions in the same basis set 
as the constant, s-wave OP component can be written 
as fj(k) = J2n,m=-oo A l nm cos[a(nk x + mky)}, provided 
that we choose the A l nm to satisfy the A\ (s-wave) sym- 
metry, A l nm = A l _ nm — A % n _ m — A % _ n _ m — A l mn , of the 
C^v character table (Table I). As for the angular momen- 
tum representation, it is then possible to write each OP 
eigenfunction in terms of the nearest pairing separations 
consistent with the representation group symmetry, mul- 
tiplied by one of the allowed generalized s-wave lattice 
eigenfunctions. 



OP Eigenfunction Crystal Representation 
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TABLE III. Lattice representations of the singlet super- 
conducting OP eigenfunctions for a tetragonal crystal. The 
functions fi(k) = X)" m= _ TO A-nm cos[a(nk x + mk v )], where 
A i — A* — A i —A* — A i 



Those OP components (or basis functions) within the 
same basis set, which have the same eigenvalues in the 
group character table, are compatible, and can in princi- 
ple mix freely, the actual mixing amount depending upon 
the particular pairing interaction form. OP components 
arising from different basis sets, with one or more differ- 
ent eigenvalues in the group character table are incompat- 
ible, and do not ordinarily mix. Each OP eigenfunction 
is a linear combination of all of the compatible OP basis 
functions within the same basis set. For tetragonal crys- 
tals, the OP eigenfunctions are listed in Tables II and 
III. 

In Table IV, we list the spin-singlet character table 
and group theoretic notation for the C 2v group opera- 
tions for the orthorhombic crystal structure appropriate 
for YBCO and Y124. || We also listed the angular mo- 
mentum (l) and lattice (n, m) representations of the two 



OP eigenfunctions, |s + d x 2_ y 2) and \d xy + g X y(x 2 -y 2 )) ■ 
Each of these OP eigenfunctions contains two basis sets 
of the tetragonal OP eigenfunctions in Tables II or III. 
Thus, |s + d x 2- y 2) contains an arbitrary mixing of the 
tetragonal |s) and \d x 2_ y 2) OP eigenfunctions, for exam- 
ple. 
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TABLE IV. Singlet superconducting OP eigenfunctions in 
the angular momentum (I) and lattice (n, m) representations, 
their group theoretic notations, and character table for the 
orthorhombic point group Civ in the form appropriate for 
YBCO. 
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We remark that even with only a slight (about 2%) 
crystallographic difference between the a and b axes in 
YBCO, the mixing of OP components that would have 
been in different representations (or basis sets) of the 
tetragonal group C± v , can be very substantial. For ex- 
ample, there could be a very large mixing of s-wave 
and d x 2_ y 2-w&ve OP components in YBCO. Especially as 
the penetration depth within the a6-planes measured in 
YBCO is anisotropic by about a factor of 2, and the c- 
axis Josephson tunneling between untwinned YBCO and 
Pb reproducibly gave values of 1.5 mV, || about 25% 
of the expected result of Ambegaokar-Baratoff , Q , one 
expects the isotropic £ = OP component to be at least 
25% of the maximum OP amplitude at low T. Of course, 
the penetration depth anisotropy is easiest to understand 
by accounting for the quasi-one-dimcnsional Fermi sur- 
face of the CuO chains. 



GT OP E a dl a d2 C 2 

Eigenfunction (a a ) (05) 

Ai \s + d xy ) +1 +1 +1 +1 

7 d + \/2X^=i{»« x 
x cos[2n(0 k - 7r/4)]} 

n,m En,m=o{l"™ + ''™] X 

x cos(nk x a) cos(mk y a) 

"h [Cnm ~h Crnn\ X 

x sm(nk x a) sm(mk y a)j 

^4.2 \d x 2—y2 ~\~ Qxy(x 2 ~y 2 ) 

) +1 -1 -1 +1 

x sin[2n(</) k - 7r/4)]} 

n,m En,m=o{[''»™"'' m »l X 

x cos(nk x a) cos(mk y a) 

"h [Cnm ^mn] X 

x sin(nk x a) sin(mfcj,a)} 



In Table V, we list the spin-singlet character table and 
group theoretic notation for the C 2v orthorhombic group 
operations appropriate for BSCCO and LSCO. We also 
listed the angular momentum (£) and lattice (n, m) rep- 
resentations of the two OP eigenfunctions, |s + d xy ) and 
\d x 2_ y 2 + g xy (x 2 -y 2 )}- Each of these OP eigenfunctions 
contains two basis sets of the tetragonal OP eigenfunc- 
tions in Tables II or III. Thus, \s + d xy ) contains an arbi- 
trary mixing of the tetragonal |s) and \d xy ) OP eigenfunc- 
tions, but not of the tetragonal \d x 2_ y i) OP eigenfunc- 
tion. We note that although the ac-plane is not a strict 
BSCCO crystallographic mirror plane, the 6c-plane is a 
strict crystallographic mirror plane (0^2 )• This fact is 
sufficient to separate the OP sets as if the crystal were 
fully orthorhombic. 

We remind the reader of the crucial difference be- 
tween the OPs with the same group theoretic notations 
for the orthorhombic symmetries appropriate for YBCO 
and BSCCO, respectively. Although for orthorhombic 
YBCO, OP components that have s-wave and d x 2_ y 2- 
wave symmetry in the tetragonal crystal can mix freely, 
and similarly OP components exhibiting d^y-wave and 
9xy(x 2 -j,2)-wave symmetry in the tetragonal crystal can 
mix, in BSCCO this is not the case. Instead, the OP 
components exhibiting d x 2_ y 2-wave and (^(a^^-wave 
symmetry in the tetragonal crystal can mix, as can the 
OP components having s-wave and d xy -wewe symmetry 
in the tetragonal crystal representation. 

IV. GINZBURG-LANDAU FREE ENERGY 

In this section, we examine the Ginzburg-Landau free 
energy of a superconductor with one or more OP com- 
ponents. We limit our discussion to a single Cu02 layer, 
assuming equivalent OPs on all Cu02 layers. We shall 
see that there is a distinct difference between the cases 
of two compatible OP components and two incompatible 
OP components. (AGL used the terminology "mixing" 
and "non- mixing" ) . 

A. single complex OP component 



TABLE V. Singlet superconducting OP eigenfunctions in 
the angular momentum (£) and lattice (n, m) representations, 
their group theoretic notations, and character table for the 
orthorhombic point group Civ in the form appropriate for 
BSCCO. Although the adi mirror plane symmetry is only 
approximate due to the c-axis component of Q, the adi mirror 
plane symmetry is robust in most current samples. 



We first assume a single OP component. Let the pair- 
ing interaction A(k, k') = Aoy(k)y>(k'), where the ba- 
sis function <p(k) is normalized, i. e., (tp(k)\(p(\i)) = 
1, which in angular momentum representation implies 
Jo^ ^r^ 2 ^^) = 1' or m real-space representation im- 
plies (a/27r) 2 Jli a /a dk x J^ a /a dk yV 2 (k) = 1, and the OP 
eigenfunction A(k, T) = A(T)^(k), where A(T) is the 
OP amplitude. The Ginzburg-Landau free energy of a 
spatially and temporally uniform superconductor with a 
single OP component is 

F = a (T)|A(T)| 2 + /5|A(T)| 4 , (1) 



G 



where a(T) — a^(T — T c ), (3 > is a constant, and 
T c is obtained from Ao in the usual BCS approximation. 
Minimizing F, we find A = and F N = for T > T c , 
and A(T) = [a (T c - T)/(2/?)] 1/2 and F s = -a 2 (T c - 
T) 2 /(4/3) for T < T c . 



B. Two incompatible OP components 

We now consider the case of two incompatible OP 
components. For simplicity, we consider the case in 
which the pairing interaction can be written in terms 
of products of only two of the incompatible, orthonor- 
mal basis functions </?i(k) appropriate for the crystal, 
A(k, k') = Yli=i ^io^iO^fiO^)- As examples of incom- 
patible basis functions, in Hgl201, YBCO, or BSCCO, 
we could choose the simplest d^-wave and d x 2_ y 2-w&ve 
basis functions, which are v / 2sin(20k) and \/2cos(20k) 
in the ^-representation. For either Hgl201 or YBCO, 
we could also choose the simplest s-wave and d xy -w&ve 
functions (1 and \[2 sin(20k) in the ^-representation), 
and for either Hgl201 or BSCCO, we could also choose 
the simplest s-wave and d x 2_ y 2-w&ve basis functions. Of 
course, the number of possible incompatible OP com- 
ponent pairs is limitless. The OP eigenfunction con- 
sistent with the above pairing interaction has the form 
A(k, T) = Y% =1 Aj(T)^(k), where the Aj(T) are com- 
plex constants representing the two incompatible OP 
component amplitudes. 

For different A^o, the bare transition temperatures T C i 
(obtained in the BCS model assuming only one Ai ^ 0) 
of the two incompatible components are different. The 
free energy can then be written as 



F = 



Si from S2 can then be obtained by inserting |Ai| from 
Eq. (0) into the linearized equation for IA2I, leading to 



1 c2 



T, 



<:2 



1 - vi(T cl /T c2 ) 
2a 20 f3i 



(4) 



Note that T< < T, 



C 2 if e — 2(5 > 0, as in weak coupling 
(BCS) theory. Then, solving a quadratic for |Ai| 2 and 
|A 2 

perature superconducting phase S2, 



\2\ 2 below T C 2, we obtain the solution in the low tern- 



|A 2 | 2 = 



|Ai| 



aao(l -vi) (T<-T) 



2/3 2 (l - vm) 



aw 



Td - F 



<:2 



(1-^)(T<-T) 
1 - v\V 2 



V2 = 



Q 2o(e - 25) 



4> 2 =^i± ir/2 for T < T< (S 2 ) 



(5) 



Because of the tt/2 phase difference between the two OP 
components in the S2 state, this is commonly known as 
the "Ai + iA2" state. As emphasized by AGL, these 
two (or more) incompatible OP components can only be 
simultaneously non- vanishing below a second phase tran- 
sition temperature, T^. Since no second phase transition 
appears to have been observed in any of the cuprates (at 
least, not above 1 K as yet), it is rather unlikely that 
> 1 K. In addition, the Ai + «A 2 states is nodeless 
below T< 2 . 



C. Two compatible OP components 



-e|A 1 | 2 |A 2 | 2 + (5[A 2 A; 2 + A2Af], 



(2) 



where the cti(T) = aio(T — T ci ) and the e and 5 
are constants. In weak coupling (BCS) theory, /3, > 0, 
e > 0, 5 > 0, and e - 25 > 0. Writing A, = 
|A,-| exp(iif>i), we note that the only term depending upon 
the phases ipi is the one proportional to 5, which becomes 
25|A 1 | 2 |A 2 | 2 cos[2(^ 1 - V2)]. Assuming T cl > T c2 , F is 
minimized in the first superconducting (Si) regime just 
below T c , by 



\A 1 (T)\ = [-a 1 (T)/(2fr)] 1 / 2 , 

|A 2 |=0, for T<<T<T cl (Si). 



(3) 



In the low temperature (S2) superconducting phase, both 
I Ai| ^ and |A 2 | ^ 0. dF/dtpi = is satisfied when 
sin[2(f/'i — ■02)] = 0, which minimizes F when ipi — ip 2 = 
±7r/2. The temperature which separates the regions 



Now we consider the more interesting case of two com- 
patible OP components, with elements <fii(k) for i = 1,2, 
of the same basis set. In Hgl201, for example, |s) con- 
tains the simplest s-wave and g x 2 y 2-w&ve elements (writ- 
ten as 1 and v / 2cos(40k) in the f -representation) , which 
are compatible basis functions. For YBCO, |s + d x 2_ y 2) 
contains as compatible basis functions the simplest s- 
wave and d 2 ,2_ J( 2-wave functions (1 and v^2cos(20k) in 
the I- representation), and for BSCCO, the simplest s- 
wave and d xy -w&ve functions are compatible elements of 
3 + d xy ). It is easy to generalize this procedure to in- 
clude all of the (infinite number of) components in the 
relevant basis set. 

With two compatible OP components, the in- 
teraction is generally of the form A(k, k') = 

Yli j=i V^OOAyoVjOO, where Ao is a symmetric ten- 
sor of rank 2. We could use this form of A(k, k') 
to generate the Ginzburg-Landau free energy, letting 
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A(k, T) = Ei=i A i( T )^i( k ), wh ere the A; are com- 
plex constants representing the OP component ampli- 
tudes. We would then find that the free energy would 
be given by Eq. ([|), plus the terms 7(A|A 2 + AiAJ) 
and (AiA* + A^A 2 )(^i|Ai| 2 + fi 2 \A 2 \ 2 ). || In the 
BCS approximation, the T C i arise from the diagonal ele- 
ments Xuo, and 7 and the /i^ arise from the off-diagonal 
A120 = A210 elements of Ao, respectively. One could then 
diagonalize the part of this free energy quadratic in the 
Ai by a unitary transformation, a two-dimensional ro- 
tation, as noted by AGL. |}2| However, AGL did not 
emphasize the important point that this rotation mixes 
the OP components for all T < T c . In order to clar- 
ify this point, we choose instead to first diagonalize the 

pairing interaction. To do so, we let Ao = RXqR -1 , 
where R is a standard two-dimensional matrix for ro- 
tation by an angle 9 about the 3 axis. By choosing 
tan(2#) = 2Ai2o/(Ano — A 2 2o), we obtain the diagonal- 

ized Ao, and the pairing interaction can now be written 
as A(k,k') = J2s=± l Ps(^)Xso'Ps(^'), where 



-4(T)/(4/? + )+MT)|A_| 5 



tp+ (k) = ipi (k) cos 9 + (f2 (k) sin 9 
tf- (k) = — cpi (k) sin 9 + ip2 (k) cos 9, 



(6) 



and A(k, T) = £s=± A s {T)ip s (k). We then expand the 
free energy in powers of A± , and obtain 



E 

s=± 



(as(T)|A s | 2 +/? s |A s | 4 



+e|A+| 2 |A_| 2 + 6{A 2 + A* 2 + AfA 2 _) 

+(A + A*_ + A_A* + )( P+ \A + \ 2 + ^\A_\ 2 ), (7) 

where a±{T) = a±o(T — T c ±), T c ± is obtained from 

A±o = H A no + A220 ± [(Ano - A220) 2 + 4A 2 10 ] 1 / 2 } in 
the BCS approximation. 

Now, it is evident that the two original OP components 
Ai and A 2 mix in two ways. First, and most important, 
the corresponding basis functions yi(k) mix via the lin- 
ear transformation employed to diagonalize the quadratic 
part of the free energy. Since the transformed OP eigen- 
functions are A±(T)</?±(k), both transformed OP eigen- 
functions contain contributions from both compatible ba- 
sis functions. Thus, the dominant OP A + (as well as 
the subdominant A_) could be written as a "Ai + A2" 
state. Although AGL did not emphasize it, this mixing 
does not require a second phase transition. 

Second, the subdominant OP (with the lower transi- 
tion temperature, T c _) is coupled to the dominant OP 
just below T c+ via the term ^ + |A + | 2 (A + A1 + A_A^). 
Depending upon the sign of A_ will either add or 
detract from A + , but in either case, it will be in phase 
with it. The most important part of the resulting free 
energy just below T c+ is then 



2|/x+||A_| -a+(T)/(2/?+) 



3/2 



which implies 



|A_(T)| 



2|// + | /-a+(T)\3/ 



a_(T c+ )V 2/3+ 



(8) 



(9) 



which behaves as (T c + — T) 3 / 2 . Thus, the subdominant 
compatible OP can pick up a strong temperature depen- 
dence just below T c+ , modifying the relative temperature 
dependences of Ai and A2 in the "Ai + A2" state. In 
particular, it is possible for one compatible component, 
say A2, to be substantially smaller than the other just 
below T c , but comparable to the other at low T. Again, 
this can occur without a second phase transition, a point 
that was incorrectly stated in the conclusion of AGL. 

More generally, when the (symmetric) interaction ma- 
trix Ao represents pairing involving a large (or infinite) 
number of compatible basis set elements, one diagonalizcs 
it, and rank orders the diagonalized interaction strengths 
A,io- Each eigenvalue Xuo corresponds to a bare T C i value, 
with T c = T c i. Each Xuo also corresponds to an OP 
eigenfunction which is a large (or infinite) sum of the 
compatible basis set elements. The resulting free energy 
will then contain quartic terms that mix two or more 
OPs. These terms will further mix the OPs, changing 
their effective T dependences and relative weights be- 
low T c . Thus, the most general OP eigenfunction will 
have the general form given in Tables II- V, with coeffi- 
cients (such as the a n for \s) in Table II) having some- 
what different T dependences. This implies that the k- 
dependence of the OP eigenfunction changes smoothly 
with temperature. |^3| None of this compatible mixing 
involves a second phase transition. 

V. C-AXIS TWIST JOSEPHSON JUNCTIONS 

We now consider briefly the case of a HTSC Joseph- 
son junction formed by twisting bicrystal halves an angle 
4>o about the c-axis, as pictured in Fig. 3. Let us sup- 
pose that there is only one phase transition observed in 
a given material. Certainly in the vicinity of the bulk 
transition T c , it is hard to imagine that a transition to 
a second state with a combined OP that mixes incom- 
patible components would be present, except in the case 
of the triplet OPs. However, since \(p x ,Py)} is a two- 
dimensional representation of the tetragonal group, it 
gives rise to a ll p x + ip y " (or equivalently, a ll pdi + ip<i2 ) 
state, as in the equatorial plane of the A phase of He 3 . 
In an orthorhombic crystal, the transition temperatures 
T c i,T C 2 for the two incompatible components are split, 
and a second phase transition appears at Tc2 < , below 
which the second component becomes non-vanishing. In 
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any event, the combined OP eigenfunction will be node- 
less and rather isotropic at low T. If this were true, one 
could not presume that the ARPES, penetration depth, 
Raman, and tunneling experiments were in any way prob- 
ing the superconducting OP, so that one would need to 
explain those and related experiments from a model that 
did not contain any nodes of the OP. 

Thus, we assume that we are sufficiently close to T c 
so that only one OP eigenfunction with compatible ba- 
sis elements need be considered. Let there be two crys- 
tal half-spaces, which are rotated an angle 0o with re- 
spect to each other. On the upper half-space, we let the 
OP be Au (0k — 0o/2), and on the lower half-space, it 
is A/, (0k' + <Po /2). We assume that there is sufficient 
Josephson coupling across the junction for a critical cur- 
rent to traverse it. 

As we discussed in detail elsewhere, there are two 
mechanisms by which quasiparticles can tunnel across 
the twist boundary. |p4| These are coherent and inco- 
herent tunneling, respectively. In the coherent tunneling 
process, the momenta parallel to the twist junction are 
conserved during the tunneling, whereas for incoherent 
tunneling, they are ordinarily taken to be completely ran- 
dom. We assume that the spatial average of the second 
order quasiparticle tunneling processes across the twist 
boundary can be written as 

(t(k - k')**(k' - k")) = S(k - k")[\J\ 2 S(k - k') 



+/ir 



-M], (io) 



where the overall two-dimensional <5(k — k") function in- 
sures translational invariance after averaging. The terms 



proportional to |J| and f- u 



6k') are the coher- 



ent and incoherent part of the tunneling. Ordinarily, 
one expects /inc(0k — 0k') to be a constant, as for the 
standard s-wave scattering calculation of Ambegaokar- 
Baratoff. (42) However, to allow for a small amount of 
forward scattering, or incoherent tunneling of non- s-wave 
character, we let 



oo 
£=0 



cos[lQ k - ^k')] 

2TTT ±e N 2D {0) 



(11) 



where N2d(0) — m/(2n) is the two-dimensional single- 
spin quasiparticle density of states. |24|] 

We remark that the question of coherent versus inco- 
herent tunneling is an interesting one, especially when 
one considers tunneling between two layers that are 
twisted about the c-axis with respect to each other. The 
important point is that a quasiparticle tunnels from the 
Fermi surface on one layer to the Fermi surface on the 
neighboring layer. This is true in both the normal and 
superconducting states. For the simple case of a quasi- 
two-dimensional free-particle band, the cross-section of 
the Fermi surface in the k x /k y plane is a circle centered 



about the high-symmetry T point at the center of the first 
BZ. Even for the case of layers twisted about the c-axis, 
the circular Fermi surface cross-sections are identical in 
both first BZs twisted an angle 4>o with respect to each 
other. In the coherent single particle tunneling process, 
one can in principle tunnel from one layer to the next 
from any position within the first BZ, and for such a free 
particle Fermi surface, coherent normal state quasipar- 
ticle tunneling could occur with equal probability over 
this circular cross-section, even for arbitrary 0o, &s long 
as the applied voltage V was effectively zero. 

In the cuprates, however, the Cu02 bands generally 
have the tight-binding, rather than the free-particle form. 
For BSCCO, the actual Fermi surface is rather compli- 
cated, as pictured in Fig. 4a. First of all, the primary 
Fermi surface is very similar to the tight-binding Fermi 
surface indicated by the thick curves in Fig. 4a. These 
curves were calculated using the two-dimensional tight- 
binding dispersion 



£(k) = t[cos(k x a) + cos^a)] 

—t' cos(k x a) cos(k y a) — /i, 



(12) 



where the primary Fermi surface is obtain by setting 
£(kp) = 0. In Fig. 4a, we used the values t = 2t' = 0.8/x, 
which gives a primary Fermi surface similar to that calcu- 
lated for BSCCO and observed in ARPES experiments. 

ppr 



In addition, there are the secondary Fermi sur- 
faces, which arise from the periodic lattice distortion 
Q = (0,0.212,1) in terms of the reciprocal lattice vec- 
tors. ]34| These are given by £(kp ± Q) = 0, which are 
indicated by the thin solid and dotted curves in Fig. 
4a. Note that for this two-dimensional dispersion, the 
c-component of Q is irrelevant to the Fermi surface struc- 
ture. For c-axis tunneling between adjacent layers that lie 
directly on top of one another, one could in principle still 
have coherent tunneling at any position on the multiple 
Fermi surfaces within the first BZ. When adjacent layers 
are twisted an angle 0o about the c-axis with respect to 
each other, as pictured for 0o — ir/4 in Fig. 4b, how- 
ever, the Fermi surfaces on adjacent layers intersect only 
at a finite set of points. This intersection has measure 
zero relative to the entire length of the two-dimensional 
Fermi surface. We thus expect the amount of coherent 
tunneling between layers twisted by a sizable angle 0o to 
be vanishingly small. Otherwise, preserving the wavevec- 
tor during the coherent tunneling process necessarily re- 
quires inelastic processes. Thus, except for the cases of 
0o ~ 0, 7r/2, we expect the dominant interlayer tunneling 
processes to be incoherent. 

Previously, [^4| we treated the theoretically more inter- 
esting (and more complicated) case in which the compo- 
sition of the OP components can depend upon the layer 
index, as the twist junction specifically removes transla- 
tional invariance. For the case in which the dominant 
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OP component is presumed to be d x 2_ y 2-wa,ve and the 
secondary, incompatible OP component is d xy -wave, the 
dominant OP component is suppressed by the proximity 
to the twist junction, locally reducing the suppression of 
the sub-dominant OP component in the bulk of the sam- 
ple by the dominant OP component. In addition, for all 
T < T c , this sub-dominant OP component would be lo- 
cally enhanced, although the enhancement would become 
quite weak close to T c . This would allow the local OP 
to twist, compensating for the physical twist in the junc- 
tion, and allowing a finite amount of Josephson critical 
current at all twist angles. 

We found, however, that sufficiently close to the bulk 
T c , only one OP component need be considered, except 
for the special (and only infinitesimally probable) case 
of an accidental degeneracy, whereby two incompatible 
OPs accidentally happen to have (almost) exactly the 
same bare transition temperature, and combine to give 
a nodeless Ai + iA 2 state. 0| Otherwise, assuming the 
bare T c b of the secondary d xy OP component is suffi- 
ciently small that T< B , the suppressed secondary tran- 
sition in the bulk, is unobservable experimentally, the 
Josephson coupling strength determines the T at which 
the OP twisting could occur. The weaker the coupling, 
the lower the T at which the OP can compensate for the 
junction twist by OP twisting, and the experiment be- 
comes definitive in the vicinity of T c . Such effects were 
shown in Fig. 4b of our earlier paper. j24j] In that case, we 
showed explicitly that for strong (rjd = r\ d = 1) <i-wave 
Josephson coupling across the twist junction (n d ) and 
non-junction (j]d) layers, respectively, a small but finite 
Ic{<f>o = 45°) was obtained at t = T/T c a = 0.7. However, 
reducing r\ d to 0.01 dramatically reduced I c {4>o = 45°) at 
t = 0.7. Moreover, as the experiments of Li et al. now in- 
dicate that T) = rj' << 1 in BSCCO, we have recalculated 
I c (<j>o) for the cases Vd = Vd = 0.1,0.001, and plotted the 
results for t — 0.5,0.9,0.99 in Fig. 5. The remaining 
parameters are the same as before, and are listed in the 
caption to Fig. 5. It is seen that for T c b/T c a = 0.2 
(about as high as one could possibly hope to have, in or- 
der for quasi-nodes to still be present at 1 K) , both sets of 
curves show the suppression of the dominant OP compo- 
nent near to T c , with 7 C (0 O )/I C (O) lying below | cos(20 o )| 
except for O ss 45°. They also show that I c (0o = 45°) 
decreases with increasing T/T c a and decreasing rjd = r/ d . 
Clearly, twisting of the purported ci- wave OP cannot ex- 
plain the experiment of Li et al, in which I c {4>o) wa s 
found to be independent of </>o . 

We therefore now consider the case in which there is 
only one OP, but that it may contain a complex mix- 
ture of compatible basis set functions. In the Ginzburg- 
Landau regime near to T c , we found that the Josephson 
critical current across the twist junction could be written 
in terms of coherent and incoherent parts, I^ oh ((po, T) + 



I™ h (<f> ,T) 2eiV 2D (0)6 (r)|J| 2 Im/ d0 k 

Jo 



2tt 



xAtffok- o /2)A£(0 k + 0o/2) 



and 



4 nc (0o,r)^e7V 2 2 D (O)ao(T)Im/ dfa dfa 

Jo Jo 

xA C /(0k ~ 0o/2)/inc('/ , k - 0k') 

xA^ k '+<V2), (13) 

where a (T) = vr/(4T) and6 (T) = 7C(3)/(8vr 2 T 2 ). Near 
T c , I c (<po,T) is proportional to (T c — T). 

We then have to evaluate I c {<t>o) for the various OP 
possibilities for Hgl201, YBCO, and BSCCO. To do so, 
it is convenient to define 



V ™ = e7rN 2D (0)a (T)/r ±0 

r ±0 , 26 (r)|J|V 



m 



a (T) 



(14) 



Clearly, rf nc is one-half the amplitude of the incoherent 
scattering contribution to the critical current for the s- 
wave OP component, and fji is a measure of the coherent 
and incoherent contributions to I c from the £-wave OP 
components, relative to that of 77 nc . Note that 772 oc r\ d 
in Fig. 5. 



GT 



OP 



Wo)/vh n ' 



Ai \s) 

B\ \d x 2_ y 2) 

B2 \dxy) 

M \9xy(x 2 -y 2 )) 



a l + J27=0 ijindl COS(47J0 O ) 

J2n=i mn-2b 2 n cos[(4n - 2)00] 
E~=i f]in-2cl cos[(4n - 2)0 O ] 
Hn=i V4ndl cos(4n0 o ) 



TABLE VI. Superconducting critical current across junc- 
tions twisted by <f>o about the c-axis of tetragonal crystals 
such as Hgl201 (from Table II). 



P c nc ((j) ,T), 1 24 1 where 
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In Table VI, we list the resulting I c (<j) )/r)™ c f° r * ne 
singlet OPs in tetragonal materials such as Hgl201. In 
Table VII, I c {4'o)l r to c is presented for the orthorhombic 
materials BSCCO and YBCO. We note, however, that in 
this case the actual forms of the OP eigenfunctions in- 
volved in the respective states of those two orthorhombic 
materials are different, as they are given in Tables IV and 
V. From Table VI, the generalized d-wave states (with 
B group theoretic symmetry) are indistinguishable from 
each other. The leading 1 = 2 components of these gen- 
eralized d-wave states give rise to I c (<t>o — 45°) = 0. The 
leading £ — 4 contribution to the generalized g X y(x 2 -y 2 y 
wave state likewise give rise to I c = at <fio = ±22.5° and 
±67.5°. For completeness, the leading I = 1 triplet OP 
components give rise to I c {4>o — 90°) = 0. These simple 
leading £ cases give rise to I c (4>o) |cos(&fo)|- These 
simple cases are pictured in Fig. 6. 

Including compatible components of higher £ values 
will not change these qualitative results. In addition, if 
the tunneling across the twist junction is entirely inco- 
herent, with s-wave incoherent tunneling only (t±i — > oo 
for £ ^ 0), then these states will give a vanishing critical 
current for all twist angles <fro (even for the untwisted case 
4>o = 0), and the generalized s-wave state would result 
in a constant 7 c (0q), the constant being a measure of the 
£ — 0, pure s-wave part of the OP eigenfunction. 

On the other hand, if the relative contribution of the 
pure s-wave form to the entire OP eigenfunction were 
very small (i. e., nearly vanishing), then a more com- 
plicated scenario could develop. For instance, if ei- 
ther higher ^-incoherent or coherent intertwist tunneling 
were present and substantial, and if the contribution to 
the intertwist tunneling arising from the g x ^ y 2 compo- 
nent were larger than the contribution arising from the 
s component, and the £ > 8- wave components of the 
OP vanished, then there could be angles 4>q satisfying 
4>q — j cos _1 [— do(l + fjo)/Vi a i] a t which I c (<f>o) = 0- In- 
cluding higher £ components of the generalized s-wave 
OP could cause the actual 4>q values to deviate from this 
simple formula, however, as pictured in Figs. 7a, 7b, and 
7c. Such a scenario is mainly expected in the case of a 
nearly vanishing £ = OP component, since the tunnel- 
ing at a finite twist angle is most likely incoherent for 
tight-binding Fermi surfaces. 
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TABLE VII. Superconducting critical current across c-axis 
junctions twisted by <f>o, for orthorhombic YBCO or BSCCO, 
from Tables IV and V. See text. 



For an orthorhombic triplet superconductor, a twist 
angle of ±90° should give rise to a vanishing critical cur- 
rent, as pictured in Fig. 6. This is clearly contradicted 
by the present experiments on BSCCO of Li et at, which 
show that a twist angle of 90° is indistinguishable from 
an untwisted junction. |^],[| Thus, it is tempting to use 
this result to rule out the triplet OPs. However, as noted 
above, the near degeneracy of the p^i-wave and pd 2 -w&ve 
OPs in BSCCO could allow the dominant OP to twist, 
giving rise to a very small, but nearly isotropic I c ((f>o). 
Thus, we cannot rule out the triplet states without some 
additional information. As mentioned above, however, 
such triplet states are expected to be nodeless for most 
of the regime T < T c . 

Now, the more interesting cases are the singlet spin 
states. There are two possibilities, states with group 
theoretic notation A\ and A?,. In Figs. 7a, 7b, and 
7c, we have pictured some examples of the various pos- 
sibilities. In the A 2 case, there is no s-wave compatible 
component, so there are quite generally twist angles </>5 at 
which the critical current vanishes. For the pure d xy -wa,ve 
OP for YBCO, or a pure d^^-wave OP for BSCCO, 
I c (c/)q = ±45°) = 0, as illustrated in Fig. 6. Similarly, for 
the OP case with GT notation A 2 , of a pure g X y(x 2 -y 2 )~ 
wave OP, I C {4>1) = at = ±22.5°, ±67.5°, as shown 
in Fig. 6. More generally, we then expect that there 
will be twist angles between these two limiting cases at 
which the c-axis critical current will vanish, as illustrated 
in Figs. 7a, 7b, and 7c. For example, if only b\ and b 2 
are non- vanishing, then the critical current vanishes at 

05 = 2 cos M L — >—), (15) 

where z = (?74 /?72 ) (^2 ) 2 ■ This expression would be 
modified somewhat by including £ = 6 and higher order 
components, but the results would still show a strong 
anisotropy of the critical current, with some angles <pQ at 
which I c {4>o) = 0- Of course, if the intertwist tunneling 
were only incoherent, but contained a small amount of 
d-wave incoherent tunneling (t± 2 < oo), but no higher- 
£-wave incoherent tunneling, then 4>q = ±45°. 

As indicated in Fig. 7a, the Ai state, which contains 
the s-wave component, can give the least anisotropic 
(and non- vanishing) I c {4>o) behavior. However, if the 
s-wave component is very small (Fig. 7b) or vanishes 
(Fig. 7c), then a highly anisotropic I c ((po) can occur. 
In YBCO, most workers apparently believe that the OP 
is \s + GJj.2_.y2), which contains d x 2_ y 2-wave and s-wave 
OP components. Although the experiment has not yet 
been performed in YBCO, this simple analysis could give 
a measure of the relative mixing of s-wave and d x i_ y i- 
wave OP components in YBCO, as a function of T. 

In addition, this experiment would allow for a mea- 
surement of the relative weight of coherent to incoherent 
tunneling, both of which are thought to be present in 
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overdoped YBCO, at least. Since coherent tunneling is 
expected to occur only for 4>q = 0, 7r/2, measuring I c {4>o) 
for <fio near to these values could give valuable informa- 
tion in this regard. For cf>Q values away from 0,7r/2, one 
could measure the relative importance of the s-wave and 
d-w&ve parts of the incoherent tunneling, along with the 
relative weight of each OP component as a function of 
T. In particular, this would give a measure of the lower 
limit of the amplitude of the s-wave component of the 
OP. To date, no such lower limit has been placed in the 
vicinity of T c for YBCO. 

In BSCCO, the situation is fundamentally different, as 
the d a; 2_ y 2-wave OP is the leading component of the A2 
state rather than the A\ state. Thus, if the present ex- 
periments of Li et al. stand the test of time, the lack of 
any <f>o dependence of I c across the BSCCO c-axis twist 
junctions is prima facie evidence that the dominant OP 
is \s + d xy ), which does not contain any component of the 
purported d x 2_ y 2-w&ve symmetry. Furthermore, it would 
be strong evidence that a d x 2_ y 2--w&ve OP component (as 
a part of the incompatible \d x 2_ y 2+g xy ^ x 2_ y 2^) OP) could 
only appear below a second phase transition. Especially 
if one were to imagine that the purported d x 2_ y 2-w&ve 
OP component were to be larger than the s-wave compo- 
nent, then one would require a second specific heat peak 
at least as large as that observed at T c . In addition, the 
incompatibility of the s-wave and d x 2_ y 2-viccve OP com- 
ponents would mean that the OP containing a purported 
d a .2_ y 2-wave component would be of the "A1-MA2" type, 
which would not have any nodes below the second tran- 
sition. 

We remark that the apparent absence of any (j)Q de- 
pendence of I c in the vicinity of T c does not necessarily 
mean that the OP is pure s-wave. It only proves that 
the OP is |s + d xy ) (Table V). If the intertwist tunnel- 
ing were pure s-wave incoherent tunneling, one could get 
a (^-independent I c for any |s + d xy ) basis set combina- 
tion, as long as the s-wave, I — OP component is finite. 
That is, the d xy -w&ve or g x 2 y 2-w&ve OP component could 
even be much larger than the s-wave OP component, giv- 
ing rise to nodes in the OP, and one could still observe 
a (/^-independent I c . However, symmetry still requires 
that a (^-independent / c (<^o) precludes the simultaneous 
presence of any d x 2_ y 2-w&ve OP component. An example 
of an OP that would be roughly consistent with the k- 
dependence of the quasiparticle gap observed in ARPES 
experiments (and also consistent with the c-axis twist ex- 
periments of Li et al.) would be the so-called "extended 
s-wave" OP, proportional to |cos(20k)|, which is a spe- 
cific example of |s) listed in Table II. This OP has nodes, 
but no da;2_j,2-wave component. 

We remark that it is really possible to do these exper- 
iments as a function of T, measuring I c (T,<fto). ||] Since 
the lack of (j>o dependence observed just below T c persists 
to low T values in all samples for which it can be mea- 
sured both in the single crystal and across the twist junc- 



tion, then it most likely that the OP is entirely \s + d xy ), 
without any d x 2_ y 2-w&ve component at any measurable 
temperature. Such a scenario would be consistent with 
the Pb/BSCCO c-axis Josephson tunneling results. Q| 
Although those authors observed very small I c R n values, 
which could by themselves be interpreted as evidence for 
a very small s-wave component to the OP, other possible, 
more mundane explanations of such small values are easy 
to imagine, even if the superconducting OP were essen- 
tially pure s-wave. [Q For example, the same low I c R n 
values were observed in NCCO, || which suggests that 
problems of a materials nature, such as an impedance 
mismatch in c-axis tunneling between cuprates and Pb, 
might be responsible for reducing I c R n . 

VI. IMPROVED TRICRYSTAL EXPERIMENT 
PROPOSAL 

We note that the tricrystal experiments of Kirtley, 
Tsuei, and collaborators have not been confirmed in a 
second laboratory. This could be due to the fact that 
it is very expensive to prepare the SrTi03 tricrystal 
substrates upon which the cuprate films were grown. 
Such depositions led to many complications at the grain 
boundaries, as have been evidenced by detailed transmis- 
sion electron microscopy (TEM) in the case of YBCO 
junctions. jl5| However, in the cases of other cuprates 
there have not been any TEM measurements available to 
the general public. Presumably, this is due to the large 
substrate costs involved. 

We thus propose a new procedure, which does not em- 
ploy any expensive tricrystal substrate, but only requires 
an inexpensive insulator to support the tricrystal. The 
improved geometry for the "tricrystal" experiment is ac- 
tually constructed out of a tetracrystal ring, as pictured 
in Fig. 8. In this setup, a single crystal of BSCCO is 
cleaved twice, giving three pieces, which we illustrated 
as dark, intermediate, and light. The intermediate piece 
is cut normal to the layers into two pieces of identical 
thickness. The dark crystal is placed upon the support. 
A section of support substrate of thickness equal to that 
of the dark crystal is placed aside it as a support for the 
rest of the ring. Then, the two pieces of the interme- 
diate crystal are placed across the dark crystal and the 
support, forming angles (^12 and (^23, as pictured. The 
two pieces of the intermediate crystal do not touch each 
other. Instead, the connection between them is made 
by placing the light crystal on top of both intermedi- 
ate crystals, forming a straight (0° or 180°) angle with 
one of them. Then, the entire triangular ring is fused 
together just below the melting point, and mounted on 
the insulating support. The experiment to test for the 
flux trapped within the tricrystal can now be performed, 
whether by measuring the flux within the tetracrystal 
with a SQUID microscope, or by applying leads, and 
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measuring the tetracrystal SQUID characteristics. In 
scanning with a SQUID microscope, one likely needs to 
fill the central region of the tetracrystal with an insulat- 
ing material such as epoxy, and to make the triangular 
central region small enough to enhance the microscope 
sensitivity. 

Let us label the crystals (1), (2), (3), and (4), clock- 
wise beginning with the intermediately-shaded crystal in 
the lower right-hand part of Fig. 8. Then, the angles 
of the junctions as one goes around the tetracrystal are 
012, 023, 034 = 0° or 180°, and 04i = ±03i- For each 
of the four junctions in this tetracrystal, the critical cur- 
rent can be calculated for the OP eigenfunction appro- 
priate for the BSCCO crystal symmetry, taking into ac- 
count the Josephson coupling strengths. If the s-wave 
OP component dominates, as in the experiment of Li et 
at, I c would be predicted to have the same sign across 
each of the junctions, so that the number of 7r-junctions 
would be 0. However, if the d x 2_ y 2 OP component domi- 
nates, then the critical currents behave as I c \2 cos(20i 2 ), 
/ c23 cos(20 23 ), ±Io34, and ±I c41 cos[2(0i 2 +02a)], where 
the ± are consistent, and relate to whether it is a 0° 
or 180° junction. For choices of the angles such that 
cos(20i 2 ), cos(20 23 ) and cos[2(0i 2 + 23 )] are all < 0, 
then it is easy to show for the 16 possible configurations 
that one always has an odd number of 7r-junctions. Thus, 
if the d x 2_ y 2--w&ve OP component dominated the c-axis 
Josephson tunneling, as required for a bulk d-wave su- 
perconductor, then one ought to have a half-integral flux 
quantum trapped in this ring at low temperatures. More 
complicated scenarios such as those pictured in Figs. 7a, 
7b, and 7c can be investigated by varying the tricrystal 
junction angles. Note that for each case in these figures 
for which / c (0o) = 0, the relevant I c changes sign at 0q. 
When only one junction is studied, the choice of phase 
factor across the junction is arbitrary, and that which 
minimizes the free energy of the junction changes I c to 
\I C \ when I c < 0. But, for a tricrystal (or tetracrystal), 
the relative phases across the junctions are fixed, due to 
the loop. 

The free energy of a ring with n Josephson junctions 
is given by 

n 

F = LI 2 /2 - ($ /2tt) £ 7 cM+ i cos(A^. M+1 ), (16) 
where 

n 

&<Pi,i+i = -2tt$/$ 

8=1 

LI 2 /2 = ($ - $ cx t) 2 /2i. (17) 

These were only written down for Josephson junctions in 
the afe-plane, as in the tricrystal experiments of Kirtley 
and Tsuei. |q j3^ , |47| However, they also apply to the case 



of a ring of c-axis Josephson junctions. One has only to 
be careful that the self-inductance L is sufficiently large 
that (3 = LI c /Qq >> 1, where I c is the minimum value of 
the critical currents across the junctions. |^6| For the ab- 
plane junctions, although J c is large, the junction area is 
small enough that I c ~ 2 mA might not be large enough 
relative to L. |s) Although J c for the c-axis junctions is 
much smaller than for the adjunctions, the much larger 
junction areas make I c ~ 20 — 200 mA, typically larger 
than for the afr-plane junctions. || Since L is propor- 
tional to the area inside the ring, by making the ring 
out of large single crystals, it should be possible to con- 
struct rings with f3 » 1. However, in order not to lose 
the SQUID sensitivity, the central area inside the ring 
should be nearly as small as for the afr-plane tricrystal 
experiments, with a rather flat insulating support for the 
SQUID microscope. 

We remark that in the previous (a6-plane) tricrys- 
tal experiments, the dependence of the critical current 
P L ? upon the misalignment angles 9i and 8j was as- 
sumed to be proportional to either cos(26*j) cos(2#j) or 
cos[2((9j + 6j)], for coherent or incoherent a6-plane tun- 
neling, respectively. However, all published experiments 
on YBCO showed an exponential dependence upon the 
misalignment angle 9 = 9i + 9j. |4^-|5ll| In the range of 
misalignment angles used in those afe-plane tricrystal ex- 
periments, both of the above theoretical predictions are 
inaccurate by several orders of magnitude. || 

Gurevich and Pashitskii showed recently that the expo- 
nential dependence on misalignment angle observed ex- 
perimentally in such YBCO afr-plane junctions is com- 
pletely unrelated to any purported OP symmetry, as 
they were able to fit the data equally well with s-wave 
and d x 2_ y 2-w&ve OP models, assuming the presence of 
insulating dislocation cores and OP suppression at the 
grain boundaries. [^2| Such grain boundary suppression 
is likely to be accompanied by oxygen stoichiometry vari- 
ation and hence local magnetic moments, which could 
be responsible for the 7r-junctions apparently observed 
in those experiments. |5l|] Such problems are evidently 
not present in the c-axis twist experiments, however. |3J 
Thus, our proposed new type of tricrystal (or tetracrys- 
tal) experiment is several orders of magnitude more reli- 
able than is the a6-plane tricrystal experiment. 

VII. CONCLUSIONS 

We have presented tables of the allowed OP eigenfunc- 
tions for tetragonal cuprates such as Hgl201 and NCCO, 
and for the orthorhombic cuprates YBCO and BSCCO. 
For tetragonal crystals, the possible spin singlet OPs are 
\s), \d x 2_ y 2), \d xy ), and \g xy (x 2 -y 2 )), tne eigenfunctions 
of which are listed in Tables II and III. For YBCO, the 
spin singlet OPs are \s + d x 2_ y 2~) and \d xy + g xy ( x 2^ y 2^) , 
with eigenfunctions listed in Table IV. For BSCCO, the 
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spin singlet OPs are \s + d xy ) and \d x 2_ y 2 + g X y(x 2 -y 2 )) i 
with eigenfunctions listed in Table V. 

We distinguished between compatible OP components, 
which are elements of the same OP basis set, and incom- 
patible OP components which are elements of different 
OP basis sets. Specifically, the s-wave and <ij.2_ y 2-wave 
OP components are compatible in YBCO, but incompat- 
ible in BSCCO. Compatible OP components have ampli- 
tudes which are in phase each other, and do not require 
the occurrence of a second phase transition to exist below 
T c . Neither do any of them require a second transition 
to become large at low T. Incompatible OP components, 
on the other hand, do not mix just below T c . Whichever 
has the higher bare T c value dominates just below T c , 
and the other can only become non-vanishing below a 
second phase transition. These incompatible OP compo- 
nents are then tt/2 out of phase with each other below 
the second transition. 

In addition, we evaluated the dependence of the o 
axis critical current upon the twist angle </>o for each of 
the compatible OP eigenfunctions appropriate for each 
of these crystal symmetries. We have discussed these re- 
sults with regard to the very recent results of Li et al. 
on c-axis twist junctions of BSCCO. ^|J3) These experi- 
ments do not show any (f>o dependence of the c-axis crit- 
ical current density just below T c . As a minimum, these 
experiments indicate that the OP is \s + d xy ). 

If the intertwist tunneling were predominantly coher- 
ent, this would further imply that the OP was very nearly 
isotropic, having an £ — s-wave form, with any addi- 
tional components of the \s + d xy ) OP being too small to 
observe. However, the more likely scenario (for BSCCO, 
especially) is that the intertwist tunneling is entirely in- 
coherent, and strongly dominated by the s-wave incoher- 
ent tunneling matrix element, such that \J\' 2 = and 
t±_o/t±£ <C 1 for £ =/= 0. In this case, fji « Sio in Table 
VII, and the c-axis twist experiments only prove that the 
OP contains an £ — 0, s-wave component, and is thus 
\s + d xy ). As such, the experiments do not rule out the 
possibility of OP nodes. Such a scenario could be con- 
sistent with the quasiparticle gap observed in ARPES, 
|p8| , ^9| with the linear, low-T in-plane magnetic pene- 
tration depth measurements, [pof and with the Raman 
scattering. [p2[ 

Assuming the OP is |s + d xy ), if the ARPES mea- 
surements were indeed giving information related to the 
superconducting OP, indicating nodes (or near-nodes) at 
0k = ±45°, then the OP would have to be similar to the 
"extended s-wave" form, |cos(20k)|. This would mean, 
of course, that the sign of the OP was constant. It would 
also mean that the amount of the d xy -w&ve OP compo- 
nent relative to the magnitude of the entire OP would be 
small. However, if the OP did not actually change sign, 
then it would be difficult to explain the impurity depen- 
dence of the in-plane penetration depth. In addition, the 
linear, low-T dependence of the out-of-plane penetration 



depth, A C (T) would be difficult to explain if the intrinsic 
interlayer tunneling were incoherent. |20| Hence, some of 
these experiments might require new theoretical expla- 
nations. 

For either predominantly coherent of incoherent inter- 
twist tunneling, the c-axis twist experiments of Li et al. 
still demonstrate the absence of any purported d x i_ y i- 
wave OP component near to T c . Since s-wave and d x i_ y i 
OP components are incompatible, any d x 2_ y 2-wave com- 
ponent could only occur below a second phase transition. 
For its amplitude to be large at low T c , the magnitude of 
the specific heat anomaly required at the second phase 
transition would have to be comparable in magnitude to 
the anomaly at T c , which can be ruled out. Thus, while 
these experiments do not preclude a small d x 2_ y 2-w&ve 
component at low T, they are inconsistent with a small 
one near to T c , and also with a large one at low T. As 
such, these new phase-sensitive c-axis twist experiments 
are in direct conflict with the recent tricrystal experi- 
ments of Tsuei and Kirtley involving BSCCO. 

It remains to be seen if the same results can be ob- 
tained in YBCO, NCCO, and Hgl201. In YBCO, since 
the s-wave and c? a ,2„ y 2-wave OP components compatible, 
one can not be as clear about the possible mixing. Never- 
theless, this experiment allows us to place limits upon the 
purported mixing of s-wave and d x 2_ y 2-w&ve OP compo- 
nents as a function of T. In no other phase-sensitive 
experiment has this been yet possible in the vicinity of 

Finally, we proposed that this technique can be modi- 
fied to prepare a tetracrystal, which can have the desired 
characteristics of the tricrystal experiments of Kirtley, 
Tsuei, and collaborators. However, in this case, the qual- 
ity of the junctions would presumably be improved by 
many orders of magnitude. In addition, they do not re- 
quire any expensive substrates, so it should be relatively 
easy and inexpensive to reproduce these experiments in 
other laboratories. 
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FIG. 4. (a) Tight-binding primary (thick solid curves) and 
secondary (thin solid and dashed) Fermi surfaces of BSCCO. 
The high symmetry points V, Y, X, and M in the first (ap- 
proximately tetragonal) Brillouin zone (BZ) are indicated, 
and the exact (r — Y) and approximate (r — X) mirror-plane 
symmetry diagonals are represented by the thin straight lines, 
(b) Two first BZs as in (a) rotated 45° with respect to each 
other. 

FIG. 5. Plot of I c (4>o) I lc(0) for the case considered in Ref. 
(24) of a dominant d x 2_ y 2-wa,ve and subdominant d^-wave 
OP, the relative amounts varying with layer index away from 
the twist junction. The parameters for these curves are 
T cB /T cA = 0.2, T< B /T cA = 0.1304, e/6(3 A = 0.5, S/6(3 A = 0.1, 
and the curves for 77 = 7/ = 0.1 and 77 = 7/ = 0.001 are indi- 
cated. Curves for t = T/T cA = 0.99, 0.9, 0.5 are presented. 



FIG. 6. Plot of J c J (0 o )/J c s versus 20o/7r for the simplest 
cases of pure s-wave (s), either d-wave (d), either p-wave (p), 
or either g-wave (g) OP contributions. 



FIG. 7. Theoretical J c J (0o)/J c S ver- 

sus 2(j>o/n for mixed OP contributions, (a) Ai state with 
[A + Bcos(2<t) ) + C*cos(40„)]/(A + B + C), and B/A, C/A 
ratios given, (b) Same as in (a), but with small A values, (c) 
A 2 state with [B cos(2^>o) + C cos(4</> )]/(.B + C), and B/C 
ratios given. 

FIG. 8. Proposed configuration of a c-axis version of the 
tricrystal ring experiment. Dark crystal: bottom. Light crys- 
tal: top. Intermediate shading: equal thickness crystals. Ar- 
rows indicate the direction of a given single crystal axis. 



FIG. 1. Two-dimensional representations of a single Cu02 
plane, (a) Idealized positions in a tetragonal crystal of the Cu 
(solid) and O (open circles) ion positions, (b) Group opera- 
tions of a tetragonal crystal, (c) Group operations of YBCO. 
(d) Group operations of BSCCO 



FIG. 2. I < 4 Even angular momentum OP basis functions 
appropriate for Hgl201, BSCCO and YBCO. 



FIG. 3. Illustration of a c-axis twist junction with twist 
angle <f> . 
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